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ABSTRACT 

We present observations of 8 Galactic Bulge microlensing events taken with the 1.0m 
JKT on La Palma during 2000 June and July. The JKT observing schedule was opti- 
mized using a prioritizing algorithm to automatically update the target list. For most 
of these events we have sampled the lightcurves at times where no information was 
available from the OGLE alert team. We assume a point-source point-lens (PSPL) 
model and perform a maximum likelihood fit to both our data and the OGLE data to 
constrain the event parameters of the fit. We then refit the data assuming a binary lens 
and proceed to calculate the probability of detecting planets with mass ratio q = 10 -3 . 
We have seen no clear signatures of planetary deviations on any of the 8 events and 
we quantify constraints on the presence of planetary companions to the lensing stars. 
For two well observed events, 2000BUL31 and 2000BUL33, our detection probabilities 
peak at ^30% and ~20% respectively for q = 10~ 3 and a ~ i?g for a Ax 2 threshold 
value of 60. 

Key words: Stars: planetary systems, extra-solar planets, microlensing - Techniques: 
photometric - 



1 INTRODUCTION 

Microlensing alters the path followed by the photons emit- 
ted by a background stellar source as they come near the 
influence of the gravitational field of a massive foreground 
object which acts as a lens. The separation of the images 
created by the lensing effect (~ 10 _3 arcsec) is too small to 
be resolved and only the combined flux is observed. The re- 
sulting lightcurve is symmetric in time with its maximum 
amplification at the time of closest approach between the 
projected position of the source on the lens plane and the 



Bascc 



lens itself. Its shape is well described by the formula: 



A(u) = 



u z + 2 



(1) 



U(u 2 +4)V2 

where A(u) is the total amplification and u is the angular 
separation of so urce and lens in units of the angular Einstein 
ring radius 9^ (Paczynski 1986). 

Even though most of the cases can be adequately de- 
scribed by this simple model, the shape of the lightcurve 
may not be symmetric and may be exhibiting significant 
deviations. These so-called anomalies of the lightcurve can 
be due to several factors and have been extensively exam- 



ined in recent literature ( Dominik 199S Wozniak . 
19971: 



ski 



scope o perated on the island of La Palma bv the Isaac Newton 
Group in the Spanish Observatorio del Roque de los Muchachos 
of the Instituto de Astrofisica de Canarias. 



1995 



Buchalter fc Kamionkowski 1997: Alcock et al 



Paczyn- 



Gaudi fe Gould 1997). Arguably the most interesting 



of these are the anomalies which can be attributed to the 
binary nature of the lens. The possibility of the secondary 
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component being an object of planetary characteristics has 
spaw ned dedicated observing campaigns to reveal their pres - 



readout, -.- photon, — sky, - total 



ence ( Albrow et al. 1998 



Rhie et al. 200C 



Bond et al. 2001) 



2 OUTLINE OF THE OBSERVATIONS 

From 2000 June 6 to July 17 we used the lm Ja- 
cobus Kapteyn Telescope (JKT) on La Palma (Longitude: 
17°52'41" West, Latitude: 28°45'40" North), Spain, to ob- 
serve a number of microlensing events alerted by the OGLE 
(Optical Gravitational Lensing Experiment) team. A CCD 
SITE2 chip was mounted on the telescope with a gain value 
of 1.95 e"/ADU and a readout noise of 11.7 e~ in TURBO 
readout mode, linear to ±1.1% for exposures to 55k ADU. 
The pixel scale was 0.33 arcsec/pixel and the observations 
were taken using an I filter at 820nm where the quantum 
efficiency of the chip is ~60%. 

Since the Galactic Bulge is a southern object, most of 
the observations were performed at high airmass (~ 2) while 
our typical seeing ranged from 1 to 2 arcsec. This however 
does not prove to be as great a hindrance as one may expect 
since photometric accura cy down to ~ 1 — 2% is still possi- 
ble [(Tsapras et al. 2001) and Fig. |l) for the brighter part 
of the lightcurves (7 < 16 mag). Exposure times were varied 
depending on the target current predicted brightness in the 
/-band so as to maximize the S/N while avoiding saturation. 
Two exposures were taken on each night for the brighter tar- 
gets to permit discrimination of cosmic ray hits. For seven of 
the events followed we obtained data at times when OGLE 
lacked any. In section 3 we discuss our observing strategy. 
Section 4 of the paper deals with the photometric analysis of 
the data. We present the analysis of the ligh tcurves in sec - 
tion 5 together with the OGLE information (Udalski 1994). 
Section 6 discusses the limits on planetary companions for 
the events. We summarize our results in section 7. 



3 OBSERVING STRATEGY 

Since our targets are only observable for about 3 hours per 
night from the JKT site during the summer, and since we 
were allocated only 2 hours per night for the observations, 
we had to fine tune the observations in order to maximize 
the number of events followed per night and also not to miss 
any events that were nearing maximum amplification. 

For this purpose, we set up a daily auto-generated web- 
page to drive the observing schedule and which the observers 
could consult at any time. The events available on each night 
were listed, along with their RAs and Decs, predicted cur- 
rent I magnitude, suggested exposure times, finder charts, 
current lightcurve fits and first estimations of the parame- 
ters. A (dimensionless) priority number was assigned to each 
event on a daily basis and the events available on each night 
were sorted and observed in order of importance. The prior- 
ity number for each event was calculated using an empirical 
formula that places more weight on short, bright, high am- 
plification events that are approaching their peaks. 

Use of a prioritizing algorithm to optimize the follow-up 
observing of microlensing events will become more essential 
from 2002, when OGLE III comes online, since the number 




Figure 1. Rms values versus the corresponding Magnitude val- 
ues of 17 measurements for all the stars found on the field of 
2000BUL34. The dominant source of noise for / > 17mag is from 
the sky. 



of microlensing events alerted is expected to rise from the 
current rate of -~100 to perhaps 300-1000 events per year. 

As automatic data reduction pipelines will soon be able 
to detect anomalous data points in real time {\ 2 of single 
lens fit > X 2 thr)' subsequent observations can be immedi- 
ately dedicated to verifying and recording the shape of the 
lightcurve anomaly. 



4 PHOTOMETRIC ANALYSIS 

The photometric analysis was perform ed using the 
IRAF/DAOPHOT package ( |Davis 199^ ) in a semi- 



automated pipeline. For further processing and lightcurve 
fitting we used separate programs developed in IDL (Inter- 
active Data Language). 

Standard pre-processing, which involves subtraction of 
the bias level and division by the master flat field, is first 
applied to the CCD frames. Once this is complete, we reg- 
ister the frames to the nearest pixel using FIGARO so that 
the star positions correspond to the same pixel areas on 
each frame. We then crop the frames to 300 x 300 pixels 
(99" x 99"), centered on the target, to reduce processing 
time since we are mainly interested in the target star and 
a certain number of nearby reference stars to determine the 
frame-to-frame magnitude corrections. The initial average 
sky value and its variance for each frame, which are needed 
to drive DAOPHOT, are determined by picking a star-free 
region and calculating the sky background photon statistics. 
The sky value is then recalculated within the routine. 

The next step involves identifying all the stars on a good 
seeing frame by using the DAOFIND task. A number (~ 
10-15) of bright, uncrowded stars are manually selected as 
point-spread function (PSF) stars which are used to create 
the PSF for the frame. The profile uses the 'penny2' analytic 
function which comprises of a Gaussian core with Lorentzian 
wings. 

We then run ALLSTAR on our frames to measure all 
the stars present in the list. The magnitudes of the PSF 
stars are used to calculate the average magnitude offset of 
each star from its weighted average over time. The average of 
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Figure 2. x 2 minimization fits to the combined data. The squares indicate JKT data and the diamonds OGLE data. 
Table 1. Fitted parameters for the events observed and reduced x 2 values. 



EVENT 


^JKT 


^OGLE 


7 OGLE 


AI 


.40 


*E(d) 


to(HJD 245+) 


x 2 

(AT-6) 


/ 


cro (mag) 


x 2 

(JV-8) 


bo 


BUL26 


21 


217 


14.04 


-0.150 


1.59 


65.60 


1705.72 


5.29 


1.939 


0.009 


1.08 


0.009 


BUL29 


5 


76 


17.07 


-0.689 


3.82 


47.43 


1698.26 


2.21 


1.000 


0.023 


1.08 


0.470 


BUL31 


29 


68 


15.49 


+0.166 


7.55 


55.79 


1743.27 


3.28 


1.738 


0.001 


0.93 


0.001 


BUL33 


26 


77 


16.99 


-0.537 


7.24 


74.26 


1730.44 


9.47 


1.875 


0.031 


1.36 


0.026 


BUL34 


17 


72 


17.73 


+0.287 


5.23 


27.04 


1710.71 


2.25 


1.441 


0.002 


1.11 


0.030 


BUL36 


18 


69 


16.06 


-0.143 


1.25 


49.03 


1711.66 


2.67 


1.628 


0.001 


1.04 


0.900 


BUL37 


16 


48 


17.62 


+0.712 


2.24 


33.08 


1713.88 


1.98 


1.227 


0.012 


1.25 


0.900 


BUL39 


13 


61 


17.11 


+0.712 


1.60 


34.87 


1719.34 


3.25 


1.679 


0.004 


1.19 


0.888 



the offsets is the zero-point value 2». This frame-dependent 
zero-point is then applied to the magnitude values of all 
measured stars on every frame. Analytically, the weighted 
average value of star nij, for each frame i, is calculated using: 

^2 mji Wji 
<m j >=^— , (2) 

i 

where the weights used are inverse variance weights Wji = 



l/a 2 ji, mji is the magnitude value of star j on frame i and 
the zero point for frame i is: 

X ( m J' 4 ~ ' m i >) w i i 
Zi = ~ ^ • (3) 

3 

Instrumental magnitudes Iji — niji — Zi rescaled in this fash- 
ion were then plotted vs time t to produce the lightcurves. 
The typical accuracy obtained by our method is illus- 
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Figure 3. Fit residuals of the combined data. The squares indicate JKT data and the diamonds OGLE data. 



trated on figure ^ which shows, for all stars found on the 
frames of the field of event 2000BUL34, the RMS scatter 
about the weighted average of 17 calibrated magnitude val- 
ues. 



5 POINT-LENS FITS TO LIGHTCURVES 

An effect that we expect to be affecting all lightcurves to 
some extent is blending. Blending is common in the pho- 
tometry of crowded fields such as the Galactic Bulge. Its in- 
fluence might lead to misinterpretation of the baseline mag- 
nitude of the source and thus inaccurate estimation of the 
true maximum amplification and timescale of the lensing 
event. Blended events may be chromatic so multiband pho- 



tome t ry can help est imate this effect (Wozniak . 
1997| ; |Vermaak 2u"0(i| ). 



Paczynski 



We performed a combined \ minimization fit to both 
our dataset and the OGLE publicly available data by op- 
timizing 6 parameters: the time of maximum amplification 
to, event timescale tg, maximum amplification Ao, baseline 
magnitude for the OGLE data -?OGLE> * ne magnitude off- 
set AJ = ^OGLE — ^JKT between the two datasets and the 
blend fraction bo = fa/ fs- fa is the blend flux of unresolved 
light sources and f 3 the flux of the unlensed source. Then the 



total observed flux at time t becomes fa Q ^ (t) — f s x A(t) + fa, 
where A(t) = (u 2 (t) + 2)/{u(t)^Ju 2 (t] )+ 4). The observed 
magnification in this case is 



A obs (t) 



fa x A{t) + fa_ A(t) + b 



(fa + fa) 



(l + 6o) 



(4) 



Our 6 parameter fits to the lightcurves of the events 
give x 2 /{N — 6) values that are generally higher than their 
expected value 1 ±(2/(iV — 6)) 1,/2 . To remedy this, we refit 
the combined dataset by introducing two additional param- 
eters that adjust the error bars, / and ao- / respresents a 
scale factor for the error bars accounting for the under- or 
over-estimation of the initial error bars. The ao term is an 
additive flux error intended to account for crowded field ef- 
fects and is most important when the source is unmagnified. 

So the error bar Si on the flux assumes the form: 



I 2 i -f 2 2 
(a + J at 



1/2 



(•») 



where the initial estimates for ao and / are and 1 respec- 
tively. By having added these two extra parameters to adjust 
the size of the error bars, we cannot use \ 2 minimization to 
optimize the fit since the fit can achieve \ 2 — by making 
the error bars infinitely large. Therefore we use a maximum 
likelihood criterion. Assuming Gaussian error distributions 
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Figure 4. x 2 minimization fits to the combined data. The squares indicate JKT data and the diamonds OGLE data. 



this is equivalent to minimizing 

JV 

X 2 +2^1n( Sl )- (6) 
i=i 

Thereby we introduce a penalty in making the error bars 
large since that increases the value of the second term in 
the equation. 

Our x 2 fits to the lightcurves of the events using 8 pa- 
rameters are summarised in Table |l] and plots of the fits are 
presented in Figures || and ^. The residuals of the fits are 
shown in | and |. OGLE data are indicated by diamonds 
and JKT data by squares. Table [I] gives for each event the 
number of JKT and OGLE data points fitted and the 8 pa- 
rameters arising from the fit. 

For 2000BUL26 we obtained data during the second half 
of the lightcurve. The scatter in the data is generally con- 
sistent with noise. No significant deviations (> 5<r) from the 
point-source point-lens (PSPL) lightcurve are seen with the 
exception of one point at ~ 14mag which deviates by +7a. 
We cannot confirm the cause of this deviation since we have 
no access to the OGLE frames but it is unlikely to be due to 
a planetary companion since further observations obtained 
by the OGLE team close to that time do not deviate from 
the expected values. 



Because at zenith distance (ZD) >45 degrees in the west 
the JKT optics show irregular jumps that usually appear 
as double images on the CCD, half of the observations for 
event 2000BUL29 had to be rejected as inappropriate for 
photometric analysis. As a consequence, this event was not 
well sampled, and we only obtained 5 data points. The fit is 
therefore mainly constrained by the OGLE data points. 

The best sampled events were 2000BUL31, whose de- 
cline we missed because of the expiry of our allocated ob- 
serving time, 2000BUL33 and 2000BUL34 where our data 
cover a significant portion of the amplification. The good 
coverage and high amplifications of these events allow us 
to place strict constraints on the fitted parameters and to 
exclude the presence of planetary conpanions in the lens- 
ing zone 0.6 < a/ Re < 1.6 (a being the planetary orbital 
radius) with high levels of confidence as discussed in section 
6. 

2000BUL37 was again covered in the decline and we 
obtained good coverage of the second half of the peak. The 
OGLE dataset lacks any points in the decline and it is the 
JKT data that help to define the shape of the lightcurve. 
2000BUL36 and 2000BUL39 were low amplification events 
selected by the priority algorithm mainly because they were 
close to maximum amplification while the remaining ongoing 
events at the time were away from their maximum amplifi- 
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Figure 5. Fit residuals of the combined data. The squares indicate JKT data and the diamonds OGLE data. 



cation values. Clearly, the information extracted from these 
last two events is not of the highest quality as their faint- 
ness and low amplifications result in poorer data points and 
a deviation should be more pronounced to be detected. All 
data is available upon request. 



6 EVENT DETECTION PROBABILITIES 

Following the 8-para meter PSPL fit, we refit the data as - 
suming a binary lens (Witt 199C; Schneider & Weiss 1986 ) 



and proceed to calculate the net detection probability (for 
a given mass ratio q) for each of the sampled events. This 
involves two additional parameters, d and q, where d is the 
projected separation between the planet and the star, and 
q the planet to star mass ratio. A similar analysi s using a 
different method has been recen tly presented in (Saudi & 
Sackett 2000| |Albrow et al. 200l[ ) 



Prior to calculating the detection probability for a 
planet of mass ratio q we set up a fine grid of planet positions 
in x,y on the lens plane and for each of these positions we 
fit the binary model to the data optimizing all parameters. 
The density of sampling in x, y has to be dense enough so 
that planetary fits are not missed. Our grid step size spacing 
was defined as y/q/4, where q is the mass ratio. This sets up 



a very fine grid for each selected mass ratio. We then make 
a Ax 2 map versus planet position by subtracting the mini- 
mum x 2 of the PSPL fit from the minimum \ 2 of the binary 
fit for each x, y. Examples of such maps are shown in fig- 
ure 0. All x, y values where the A% 2 exceeds the threshold 
value (Ax 2 ^j lr =60) are shown in black. These 'black zones' 
show us where the PSPL model gives a better fit to the data. 
The appearance of 'white zones' on the plots would signify a 
better fit has been achieved with the assumption of a planet 
at that position interacting with one of the two images of 
the source star (A^ 2 < — Ax 2 ^ nr ). Our plots show no white 
detection zones for any of the 8 events followed. Note that 
the detection zones closer to the Einstein ring of the lens 
are larger. This is because a planet in that vicinity would be 
perturbing a more highly amplified image. 
We 

calculated the detection probability for a planet/star mass 
ratio q = 10 -3 for all 8 events followed using two detection 
threshold values Ax 2 ^ nr =25,60. The results are plotted in 
figure ^[ Briefly, the detection probability of finding a planet 
at the position x, y on the lens plane for an orbital radius a 
is calculated as: 



P(det\a) 



P(det\x, y)P(x,y\a)dxdy. 



(7) 
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(d) 2000BUL33 
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(b) 2000BUL29 



(c) 2000BUL34 



(h) 2000BUL39 



(c) 2000BUL31 



(f) 2000BUL36 



(i) all 8 events 



Figure 6. Planet detection probabilities for the events observed. The curves are for a planet/star mass ratio q = 10 3 . The results 
shown are for Ax 2 threshold values of 25 (dashed) and 60 (dot-dashed). 



The first term on the right side is 



P(det\x,y) 



1 ifA x 2 >A x 2 thr 
otherwise. 



(8) 



where Ax 2 is the change in x 2 f° r a planet at x, y relative 
to the no-planet model. This term becomes significant when 
the planet at x, y lies close to one of the images of the source 
at the time of one of the data points in the lightcurve. The 
second term P(x, y\a) is obtained by randomly orienting the 
planet's assumed circular orbit of radius a, and then pro- 
jecting it onto the x,y plane of the sky. This gives a circular 
distribution centred on the lens star and rising as id /a) 2 to a 
sharp peak at d = a, outside which the probability vanishes. 
This term may be written as: 



P{x,y\a) = 



2nay/a? 



(9) 



for d 2 = x 2 + y 2 < a 2 . A slightly elliptical orbit would blur 
out the outer edge, and it's obviously possible to calculate 
this for any assumption about the eccentricity. 

The net detection probability P(det\a) is therefore the 
result of summing up the fraction of the time that a planet 



in the orbit of radius a would be located inside one of the 
'black zones'. 

For our threshold value of Ax 2 ^ nr =60, the probabil- 
ity of detection reaches a highest value of ~ 30% for 
q = 1CT 3 and a/R E ~ 1 for event 2000BUL31 and ~ 20% 
for 2000BUL33. The remaining events that were followed 
have lower (1-10%) detection probabilities. The probability 
then drops off for larger separations. As a general trend the 
probability peaks around 1 Re as expected while for events 
with low amplifications, for the mass ratios presented, we 
observe a double peaked probabilit y in agreement with th e 
simulations performed by Gaudi (Gaudi & Sackett 200C) 
for the detection efficiency. The total probability from all 
8 events combined is obtained by summing up the individ- 
ual probabilities for each event and for each value of a /Re- 
It is dominated by the contributions from 2000BUL31 and 
2000BUL 33 a nd has a maximum value of 63% for a ~ Re 
(see Fig. 6(i)) for A X 2 ^ ur =60. This maximum approaches 



100% if we use a detection threshold to 25. 
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(a) 2000BUL31 (b) 2000BUL33 



Figure 7. Ax 2 -vs- planet position for the data on events 2000BUL31 and 2000BUL33. The black detection zones show where the 
presence of a planet with a planet/star mass ratio q = 10~ 3 is ruled out by our observations. A successful detection would have been 
indicated with a white zone. The A% 2 threshold value used for these plots was set at 60. 



7 SUMMARY 

We have followed 8 microlensing events using the JKT on 
La Palma for 2hrs per night from 6 June to 17 July 2000. 
We presented fits to the combined JKT-OGLE datasets 
and recalculated the event parameters. We searched the 
data for signatures of planetary companions with mass ratio 
q — 10 -3 but have seen no indications of a planetary pres- 
ence in the datasets. Finally we calculated the planetary 
detection probabilities on all the events for a mass ratio of 
q = 10~ 3 . For events 2000BUL31 and 2000BUL33, our de- 
tection probabilities peak at ~30% and ~20% respectively 
for q = 10" 3 and a ~ R E for Ax 2 tn r= 60 - 
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